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We present experimental evidences for the existence of a semi-hard scale in light 



t>» ■ hadrons. This includes the suppression of gluon radiation that is seen in high mass hadron 
diffraction; the weak energy dependence of hadronic total cross sections; the small value of 
the Pomeron trajectory slope measured in photoproduction of J/^f; the weakness of gluon 
shadowing in nuclei; shortage of gluons in the proton revealed by an unusual behavior 
of the proton structure function in the soft limit, and the enhanced intrinsic transverse 
momentum of quarks and gluons, which considerably exceeds the inverse hadronic size. 
All these observations suggest that gluons in hadrons are located within spots of a small 
Q_i' size relative to the confinement radius. 

1. Introduction 

^ ■ There is growing theoretical and experimental support leading towards the existence 
of a non-perturbative scale smaller than the usual 1/Aqcd ~ lfm, and which is related 
to the gluonic degrees of freedom. First, an analysis of hadronic matrix elements of the 
gluonic contribution to the energy momentum tensor, using the QCD sum rules approach, 
gives a value of 0.3 fm for the radius of the corresponding form factor [1]. From the lattice 
side, numerical simulations of the gluon two point correlation function turn out a value of 
also 0.2 — 0.3 fm for the correlation length [2], and the energy of the QCD string appears 
concentrated in a tube of radius 0.3 fm in the transverse direction [3]. On the other side, it 
has been shown that the instanton radius peaks approximately at 1/3 fm [4]. Furthermore, 
high statistics data for diffractive gluon bremsstrahlung in hadronic collisions is difficult to 
explain unless gluons in the proton have transverse momenta as high as about 0.7 GeV [5]. 
This has been confirmed by studies of diffractive parton distributions, which concluded 
that they have a rather small transverse size [6]. What actually happens is that the 
smallness of the gluon clouds slows down Gribov diffusion of the gluons in transverse 
plane, and this results in a small slope of the Pomeron trajectory in hard reactions [7], in 
agreement with data. More arguments in favor of small gluonic spots coming from DIS 
can be found in [8]. Some of these results have been corroborated by recent studies of the 
spatial distribution of gluons in the transverse direction at small x [9]. Here we overview 
the available experimental evidences for the presence of a semihard scale in hadronic 
structure. 
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2. Why gluon radiation is suppressed 

If gluons in hadrons are located within small spots of radius r , they have enlarged 
transverse momenta qx ~ l/ r o- Such gluons cannot be resolved by soft interactions and 
be shaken off, which means that the bremsstrahlung cross section should be suppressed 
compared to perturbative estimates. 

However, in the case of soft inelastic collisions followed by multiparticle production the 
events with or without gluon radiation look alike. In both cases the produced particles 
build a plateau in rapidity, and then it is difficult to find any definite signature of the 
radiated gluons. 

Diffraction offers an exclusive possibility to identify gluon radiation. A high-energy 
hadron can dissociate diffractively either via excitation of the valence quark skeleton, or 
by radiating gluons. These two mechanisms are characterized by different dependence on 
the effective mass, M x , of the excitation, 

4- excitation of the quark skeleton 

x (1) 
jr diffractive gluon bremsstrahlung 

The Mx-dependence at large Mx correlates with the spin of the slowest particle produced 
in the excitation. Only a vector particle, i.e. a gluon, can provide the 1/M X dependence. 

Thus, one can single out the cross section of diffractive gluon radiation from the large 
mass tail of the M^-distribution. An analysis [10] of diffractive data shows that gluon 
radiation is amazingly weak. In order to understand that one can interpret diffraction in 
terms of the Pomeron-proton total cross section, as is shown in Fig. I. If we treat the 



da (hp — > Xp) 
dMl 




P 



Figure 1. The cross section of diffractive excitation of a proton expressed in terms of the 
total Pomeron-proton cross section. 



Pomeron like a gluonic dipole, one may expect a cross section 9/4 times larger than for 
a qq dipole. Comparing with the pion-proton cross section, say 25 mb, one arrives at the 
estimate of about 50 mb. However, data suggest quite a smaller value, about 2mb. A 
straightforward explanation for such a dramatic disagreement would be a much smaller 
size of the gluonic dipole (Pomeron) compared to the quark- ant iquark dipole (pion). Thus, 
one concludes that gluons should be located within small spots in the proton. 

Although it is not quite rigorous, one might try to estimate the diffractive radiation 
cross section perturbatively, and in this case the result exceeds data by more than an order 
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of magnitude. To reduce the cross section down to the observed value one should assume 
that the mean quark-gluon separation is as small as r = 0.3 fm [5]. With such a modified 
quark-gluon light-cone distribution function, the effective triple-Pomeron coupling has the 
form [5], 



Here we assume that 1 3> 1 — x F 3> s /s, where s ~ 1 GeV 2 . The energy dependent 
parameters <To(s) and Ro(s) controlling the shape of the universal dipole cross section [11] 
are defined in [5]. With tq = 0.3 fm, the triple-Pomeron coupling eq. (2) agrees with the 
result of the triple-Regge analysis [10] of single diffractive data. 

3. Why energy dependence of hadronic cross sections is so weak 

It is well known that hadronic cross sections rise with energy approximately as s e , where 
the exponent is quite small, ewO.l. What is the origin of this small number? We do not 
expect any small parameters in the soft regime of strong interactions. 

This problem is closely related to the topic of the previous section. In fact, the energy 
dependence is driven by gluon radiation which turns out to be suppressed. Before we saw 
a manifestation of this effect in diffraction, now in the total inelastic cross section. 

Without gluon radiation the geometric cross section of two hadrons would be constant, 
since their transverse size is Lorentz invariant, i.e. is energy independent. The phase 
space for one gluon radiation is proportional to ln(s), so multigluon radiation leads to 
powers of ln(s) in the cross section. The calculations performed in [12] confirm this. The 
hadronic cross section was found to have the following structure, 



where <t is the energy independent term related to hadronic collisions without gluon 
radiation. The second term in (3) is the contribution of gluon bremsstrahlung to the 
total cross section. Here the parameter <j\ is expected to be small due to the smallness 
of the gluonic spots. Indeed, it was found in [12] that o\ = 270^/4, where factor 
C 2.4 is related to the behavior of the dipole-proton cross section, calculated in Born 
approximation at small separations, a(rr) = Cr\ at r T — > 0. 

The energy dependence of the second term in (3) was found to be rather steep, A = 
4a s /37T = 0.17. This exponent seems to be too large compared to the experimentally 
measured e ~ 0.1. There is, however, no contradiction due to the presence of the large 
energy independent term in (3). Approximating the cross section (3) by a simple power 
dependence on energy, the effective exponent reads, 






(3) 



A 



(4) 



1 + o Q /a l (s/s )- A 



So, one should expect a growing steepness of the energy dependence for the total cross 
section. One can estimate the value of r demanding the effective exponent to be e ps 0.1 
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in the energy range of fixed target experiments, say at s ~ 1000 GeV 2 . With <jq = 40 mb 
found in [12] one gets r = 0.3 fm. 

Thus, the observed slow rise of the total hadronic cross sections provides another evi- 
dence for the existence of small gluonic spots with transverse size r ~ 0.3fm. 

One may expect a steeper energy dependence for heavier flavors. Indeed, for J/ty- 
proton scattering <t is so small, that e k A. Indeed, data for J/^> photoproduction 
from HERA [13] show that e m 0.2. One should be careful, however, interpreting the 
data within the vector dominance model [14], and remember that Eq. (3) was derived 
assuming that r is much smaller than the hadronic size, otherwise interferences should 
be included. 



4. Why diffractive cone shrinks with energy slowly 

The prediction of a shrinkage of the diffraction cone has been one of the first achieve- 
ments of the Regge theory. Indeed, data show that the elastic slope in hadronic collisions 
rises with energy as, B e i(s) = B Q + 2a'jp ln(s/s ), where a'jp ~ 0.25 GeV -2 . This is about 
four times smaller than in binary processes mediated by other Reggeons, a' m ~ 1 GeV" 2 . 
Why? 

The diffractive cone shrinkage is usually related to Gribov diffusion of gluons in the 
transverse plane. If each "step" in impact parameters, Ab 2 = r 2 ,, is small, the diffusion 
should proceed slowly. Indeed, a rather small value of the Pomeron trajectory slope was 
predicted in [12], 

, 1 dB ei a s 2 2 
2 am(s/s ) on 

This seems to be too small, in strong contradiction with value 0.25 GeV -2 known from 
data for the elastic slope. One may wonder, why the same model [12] which predicts (5) 
describes well data for elastic slope, as is demonstrated in Fig. 2? The relatively large value 
of a'jp, turns out to result from unitarity saturation. Indeed, the elastic differential cross 
section (actually the amplitude) Fourier transformed to impact parameter representation 
(see details in [12]) demonstrate unitarity saturation at small impact parameters. In spite 
of the observed rise of the total cross section with energy, there is no room for further 
growth at small impact parameters, only the amplitudes of peripheral collisions rise with 
energy. This leads to a rising with energy radius of interaction directly related to the 
elastic slope. Thus, a substantial part of the observed energy dependence of the elastic 
slope and of the effective a' e ^ = 0.25 GeV" 2 is related to saturation of the unitarity 
bound. How to disentangle the two effects? 

To get rid of unitarity corrections one can consider the interaction of a small dipole 
with a proton. For example photoproduction of a heavy quarkonium [7], or high Q 2 
electroproduction of p. Then, the elastic amplitude is too small to be affected by unitarity 
(absorptive) corrections, and the energy dependence of the slope must be solely due to 
the rise of the gluon clouds, i.e. Gribov diffusion. This expectation of [12] was nicely 
confirmed in elastic photoproduction 7 + p ^ J/^/ +p measured by the ZEUS experiment 
[13] which found a' w = 0.115 ± 0.018 GeV" 2 . The data and fit are depicted in Fig. 3 
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Figure 2. Elastic slope for pp/pp 
(bottom/upper curves) collisions, 
calculated in [12] with the small 
Pomeron slope according to (5). The 
references to the data can be found 
in [12]. 
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Figure 3. Elastic slope in exclusive pho- 
toproduction 7 + p — > J /iff + p. Data 
are from the ZEUS experiment [13]. The 
solid curve is a fit with a'jp = 0.115 ± 
0.018 GeV -2 , in good agreement with (5). 



5. Why gluon shadowing is so weak 

Although bound nucleons in a nucleus are rather well separated, Lorentz-boosted into 
the infinite momentum frame they start overlapping at small Bjorken x. This happens 
since Lorentz contraction of the gluonic clouds at small x is much weaker than at large x. 
If the clouds originated from different nucleons overlap, they may fuse reducing the gluon 
density at small x. This phenomenon is called gluon shadowing. At first glance this might 
be a considerable effect, since gluons interact stronger than quarks. However, one should 
be cautious, since a similar naive expectation of a large cross section for gluonic dipoles 
failed when was confronted with data (see Sect. 2), and both phenomena have common 
roots. 

Even if gluonic clouds overlap in the longitudinal direction, one should make sure that 
it also happens in the transverse plane. This is not obvious if the clouds are small. In a 
nucleus the mean number of gluonic spots overlapping at a given impact parameter is, 

(n) = -rl(T A ), (6) 

where (Ta) ~ Pa §-Ra is the mean nuclear thickness; pa ~ 0.16 fm -3 is the nuclear density; 
and Ra ~ 1.14 fm x A 1 ' 3 is the nuclear radius. Notice that r is the mean diameter of a 
dipole, rather than radius, which explains the factor 1/4 in (6). 
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According to (6) even for the heaviest nuclei only very few gluonic clouds have a chance 
to overlap at a given impact parameter, (n) ~ 0.3. Therefore, the longitudinal overlap 
of gluons originating from different nucleons does not lead to gluon interaction if the 
transverse overlap is so small. This may substantially weaken all collective phenomena at 
small x including color glass condensate [15]. 

Thus, another consequence of smallness of the gluonic spots is a very weak gluon shad- 
owing, predicted in [5]. Fig. 4 demonstrates results of calculations performed in [5] by 
different methods and at different scales Q 2 ~ 1/r 2 and 4GeV 2 . Even at x as low as 1CT 4 
shadowing corrections for medium heavy nuclei are only 10%. 




0.7 J —I 

1 " 4 1 " 3 1 " 2 

X 

Figure 4. Predictions for gluon shadowing 
based on the small gluonic spot structure 
[5] . The two sets of curves are calculated by 
different methods and at different scales, at 
Q 2 ps 1 jr\ (thin curves) and at Q 2 = 4 GeV 2 
(thick curves). Each set is calculated for car- 
bon, iron and led (from top to bottom). 
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Figure 5. Gluon shadowing in cal- 
cium suggested by DGLAP analysis of 
DIS data on nuclei. The solid curve is 
from an analysis [18] done in the next- 
to-leading order approximation. The 
dashed curve is the result of a leading 
order analysis [19]. 



Unfortunately, the gluon distribution function is difficult to measure, since the main 
partonometer, deep-inelastic scattering (DIS), probes directly only quarks and antiquarks. 
The gluon distribution can be accessed only via Q 2 evolution. Although previous attempts 
to single out the gluon distribution in nuclei from DIS resulted in big uncertainties, a recent 
next-to-leading order analysis [18] turned out to be quite sensitive to gluons. Gluon 
shadowing extracted from data is depicted in Fig. 5 for calcium [18]. Dashed curve 
show the results of a LO analysis [19]. We do not compare with the gluon shadowing 
suggested in [20], since it was ad hoc. The results of the DGLAP analyses agree quite 
well with predictions [5] within the theoretical and fitting uncertainties. Other gluon 
shadowing calculations missing the effect of small gluonic spots dramatically overestimate 
the magnitude of gluon shadowing. 

Suppression of gluon radiation due to smallness of gluonic spots also leads to a sub- 
stantial reduction of the effect of color glass condensate. This is demonstrated in Fig. 6 
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where py-dependent nuclear ratios are calculated for long-range (ro ~ 1/Aqcd) and short 
range (r ~ 0.3 fm) gluons [16]. Gluonic spot structure of the proton also leads to a very 
weak Cronin enhancement predicted for RHIC (confirmed) and LHC [17]. 
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Figure 6. Normalized ratios of gluon radi- 
ation cross sections in quark-nucleus over 
quark-nucleon collisions. The ratios are 
plotted as function of gluon transverse 
momentum. Calculations [16] are per- 
formed for lead in two cases of long-range 
and short range (r = 0.3 fm) gluons. In 
the latter case the magnitude of the color 
glass condensate effect is much smaller. 
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Figure 7. Logarithmic Q 2 derivative of 
F 2 (x,Q 2 ) as function of x and Q 2 [21]. 
The observed turn over and fall at small 
Q 2 indicates at shortage of gluons. 



6. Why there is a shortage of gluons at low scale 

As far as gluons are located within small spots, it is difficult to resolve them at low scale, 
Q 2 < 4/Vg. With poor resolution the proton looks like a 3-quark system containing no 
gluons. At the same time, no changes happen at higher Q 2 which resolve distances much 
smaller than the size of gluonic spots. In this regime the gluon density is rising toward 
smaller x unless related values of Q 2 become too low to resolve the spots. This explains the 
ZEUS data [21] depicted in Fig. 7 for logarithmic Q 2 derivative of the structure function. 
This derivative suddenly drops at Q 2 below few GeV 2 showing that parton distributions 
are frozen below this scale, no evolution happens due to lack of gluons. 

7. Why transverse momenta of hadron constituents are so large 

If there was only one scale in hadronic structure, the charge radius, one would expect 
the mean momenta of hadron constituents to be of the order of the inverse radius, i.e. 
(k T ) ~ Aqcd ~ 200 MeV. This is what one would observe with a poor resolution, 
insufficient for seeing the structure of constituent quarks. However, with somewhat better 
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resolution one can see gluons whose Fermi motion is much more intensive, since they are 
confined within small spots, (kx) ~ 1/Vo ~ 700 MeV. There should be manifestations of 
such an enhanced Fermi motion in reactions and observables sensitive to the primordial 
parton momentum. 

7.1. Sea-Gull effect 

The projectile quark and gluons lose (at least partially) coherence and emerge from 
the interaction area as forward jets. Theses jets having rather small transverse momenta 
hardly can be reconstructed. However, inclusively detected hadrons carrying fraction 
z of the initial quark/gluon correspondingly have same fraction of the jet transverse 
momentum. Therefore the mean transverse momentum of inclusive hadrons should rise 
as function of Feynman xf- This effect is indeed observed in data [22] and it allows to 
measure the jet transverse momentum which results from the momentum transfer gained 
in the interaction and from the primordial intrinsic momentum of the parton. The latter 
can be extracted from data and it turns out to be rather high, 0.5 — 1 GeV, the value 
anticipated in the model of gluonic spots. 

7.2. Unintegrated gluon distribution 

Direct information about transverse momenta of gluons in a hadron comes from the 
unintegrated distribution function. There are few phenomenological distributions in the 
literature. A quite popular one [23] employs the saturated shape of the dipole cross section 
fitted to HERA data for the proton structure function at small Bjorken x, 



where Ro(x) = 0.4 fm x (x/x Q ) - 1AA ; x = 0.0003. From this distribution the mean 
transverse momentum squared is (k^) = 12/Rq(x). For instance, at x — 0.01 this is 



yj (kj) ps 1 GeV. Although the distribution Eq. (7) falls off steeply at large kr, this does 

not have a big effect on its mean value. Notice that for the small kx region an energy de- 
pendent saturated shape of the dipole cross section is more appropriate [5], which results 



in a smaller value for the mean transverse momentum of gluons, J (k T ) ps 0.6 GeV [24]. 



Both estimates of the mean transverse momentum are quite large indicating at local- 
ization of gluons within small spots. 

7.3. (t l /(t t in DIS 

Of course, not only gluons get increased transverse momenta, but also valence quarks 
if they are probed with sufficient resolution, Q 2 3> 1/tq. The ratio of the longitudinal-to- 
transverse DIS cross sections, R = <Jl/&t, is sensitive to the transverse motion of valence 
quarks (at large x), since otherwise is vanishes in LO at high energy v (or high Q 2 and 
fixed x) as R = g(x)Q 2 jv 2 according to the Callan-Gross relation [25,26]. 

In order to incorporate the effect of transverse motion of the quarks, a parametrization 
R = cQ 2 /{Q 2 + d 2 ) 2 was used in [27] to fit data, with d 2 = 0.99 ± 0.39 GeV 2 . This result 
agrees with the expected d ~ 1/ro- 

Notice that data for R x , Q 2 is well explained in the dipole approach [28] with a saturated 
cross section. This is not a surprise, since according to our discussion in the previous 
section the dipole cross section bears information about enhanced gluon momenta. 
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